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In two-dimensional topological insulators, such as inverted HgTe/CdTe quantum wells, helical 
quantum spin Hall (QSH) states persist even at finite magnetic fields below a critical magnetic field 
Be, above which only quantum Hall (QH) states can be found. Using linear-response theory, we 
theoretically investigate the magneto-optical properties of inverted HgTe/CdTe quantum wells, both 
for infinite two-dimensional and finite-strip geometries, and for possible signatures of the transition 
between the QSH and QH regimes. In the absorption spectrum, several peaks arise due to non- 
equidistant Landau levels in both regimes. However, in the QSH regime, we find an additional 
absorption peak at low energies in the finite-strip geometry. This peak arises due to the presence of 
edge states in this geometry and persists for any Fermi level in the QSH regime, while in the QH 
regime the peak vanishes if the Fermi level is situated in the bulk gap. Thus, by sweeping the gate 
voltage, it is possible to experimentally distinguish between the QSH and QH regimes due to this 
signature. Moreover, we investigate the effect of spin-orbit coupling and finite temperature on this 
measurement scheme. 
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I. INTRODUCTION 

Since the first theoretical predictions of the quantum 
spin Hall (QSH) effect in grapheneii^ and in inverted 
HgTe/CdTe quantum-well (QW) structures^ topological 
insulators and topological superconductors have evolved 
into a topic of immense research interest in recent 
yearsi^^— Shortly after those proposals, the QSH state 
has first been demonstrated experimentally in inverted 
HgTe/CdTe QWs [see Fig. [T] (a)])^“— where one can 
tune the band structure by fabricating QWs with dif¬ 
ferent thicknesses dr^ Furthermore, several other two- 
dimensional (2D) systems, such as GaAs under shear 
strain^ 2D bismuthr^ or inverted InAs/GaSb/AlSb 
semiconductor QWs with type-H band alignment,— have 
been proposed theoretically to exhibit QSH states, and 
three-dimensional analogs of the QSH states have also 
been found, both theoreticalljd^ and experimentally ; 
giving rise to the concept of topological insulators : 

Topological insulators are materials that are insulat¬ 
ing in the bulk, but possess dissipationless edge or sur¬ 
face states, whose spin orientation is determined by the 
direction of the electron momentum. In analogy to the 
helicity, which describes the correlation between the spin 
and the momentum of a particle, those spin-momentum 
locked edge or surface states have been termed heli¬ 
cal. A 2D topological insulator, synonymously also re¬ 
ferred to as a QSH insulator, and its helical edge states 
are illustrated in Fig. [T] (b): Due to their helical na¬ 
ture, those topological edge states are counterpropagat- 
ing, spin-polarized states, which are protected against 
time-reversal invariant perturbations such as scattering 
by nonmagnetic impurities, and thus also of interest 
for spintronics applications— These QSH states are 





FIG. 1. (Color online) Schematic views of the (a) HgTe/CdTe 
QWs considered in this work as well as of the (b) QSH and 
(c) QH edge states and classical bulk orbits. Here, dotted 
and crossed circles denote spin-up and spin-down electrons, 
respectively, while arrows denote their respective directions 
of motion. 


in sharp contrast to quantum Hall (QH) states, which 
emerge if a magnetic field is applied and which are chi¬ 
ral in the sense that depending on the direction of the 
magnetic field they propagate in one direction only at a 
given edge—independent of spin [see Fig. [1] (c)]. 

Following the experimental demonstration of the QSH 
effect in HgTe-based QWs, much effort has been invested 
in the theoretical investigation of the properties of 2D 
topological insulators, their helical edge states, and pos¬ 
sible applications— At the heart of the QSH state are 
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FIG. 2. (Color online) Magnetic field dependence of the states 
at fc = 0 in a finite strip of width w = 200 nm compared to 
the bulk LLs. The thinner solid and dashed lines represent 
bulk LLs for s —f and s =),, respectively. The levels of the 
finite-strip geometry are displayed by thick lines. All lev¬ 
els displayed here have been calculated for band parameters 
corresponding to d = 7.0 nm, where the effective model in¬ 
troduced in Ref. @ yields a critical magnetic field Sc ~ 7.4 T. 
From Ref. [4^. 


relativistic corrections, which can—if strong enough— 
result in band inversion )^^'^^ that is, a situation where 
the normal order of the conduction and valence bands is 
inverted and which can lead to peculiar effects such as the 
formation of interface statesi^^— In HgTe/CdTe QWs, 
the band structure of the 2D system formed in the well is 
inverted if the thickness d of the HgTe QW exceeds the 
critical thickness dc ~ 6.3 nm, whereas the band struc¬ 
ture is normal for d < dc- Hence, at H = 0, QSH states 
are absent in HgTe QWs with d < dc, but can appear for 
d > de¬ 
li a magnetic field is applied to an inverted HgTe 
QW, Landau levels (LLs) and their accompanying edge 
states arise. It has been known for a long time that be¬ 
low a critical magnetic held Be the uppermost valence 
LL has electron-like character and the lowest conduc¬ 
tion LL has hole-like character in these inverted QW 
structures 42r— In this situation, counterpropagating, 
spin-polarized states also still exist. Thus, the QSH state 
persists even at hnite magnetic held B < Be, although 
these approximate QSH states are no longer protected 
by time-reversal symmetry4^“— We will call this inverted 
regime oi B < Be, where both QSH and QH states coex¬ 
ist, the QSH regime in the following. For larger magnetic 
helds B > Be, the band ordering becomes normal (see 
Fig. [2) and only QH states can be found. This normal 
regime, that is, a HgTe QW with d > dc and B > Be or 
a HgTe QW with d < dc, will be called the QH regime 
throughout the paper. 

There has been a great amount of experimental in¬ 
terest in magneto-transport and (magneto)-optical prop¬ 
erties of HgTe-based QWs^ d^d^~ d£ and other topologi¬ 


cal insulators)^— For example, magneto-transport mea¬ 
surements near the charge neutrality point in a system 
which contains electrons and holes point to so-called 
snake states, known also from other materialsplay¬ 
ing an important role in this regime— Terahertz ex¬ 
periments have revealed a giant magneto-optical Fara¬ 
day effect in HgTe thin films— Furthermore, a reso¬ 
nant photocurrent has been observed in terahertz exper¬ 
iments on HgTe QWs with critical thickness and argued 
to originate from the cyclotron resonance of a linear Dirac 
dispersion Spectroscopic measurements also point to 
a LL spectrum characteristic of a Dirac system— 

While there are theoretical studies on magneto- 
transport,—)^ anomalous galvanomagnetism,— and 
magneto-optical effects such as the Faraday and Kerr 
effects^^— in 2D topological insulators and topological 
insulator films, our goal here is to investigate magneto¬ 
optical properties of HgTe/CdTe QWs and search for 
signatures of the transition between the QSH and QH 
regimes in inverted QWs. 

The paper is organized as follows: Following the in¬ 
troduction of the model and formalism in Sec. im the 
results for a bulk HgTe QW are discussed in Sec. mil 
while Sec. IIVI contains the main focus of this work and 
is devoted to the discussion of finite strip geometries. A 
brief summary concludes the paper. 


II. MODEL AND METHODS 


A. Model system 


For a description of the HgTe/CdTe QWs situated in 
the a;y-plane and subject to a perpendicular magnetic 
held B = Bbz (with B > 0 throughout this paper), we 
use the gauge A(r) = —ByOx for the magnetic vector 
potential, which is convenient if the system investigated 
is conhned in the y-direction. Then, the QW is governed 
by the 2D, effective 4x4 Hamiltonian^)iS 


Hq =Cl4 -f AlFs 


DI4 + HF5 
¥ 
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( 1 ) 

where Px and Py are the momentum operators. A, B, 
C, T), and M. are material parameters depending on the 
QW thickness d (along the ^-direction), Ib = ^yh/e\B\ = 
y^h/eB and Pb denote the magnetic length and the Bohr 
magneton, respectively, e = |e| is the elementary charge, 
and h is Planck’s constant. The effective Hamiltonian (HD 
captures the essential physics in HgTe/CdTe QWs at 
low energies and describes the spin-polarized electron-like 
{E) and heavy hole-like (H) states jift); 1^ t)i lAj.), 
and \H j,) near the F point. 

The matrices in Eq. (HD are given by the 4x4 unity 
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matrix I4 and 


ri = 



0 


T 2 = 



rs 


(Jz 0 \ _ ( (Jg 

0 a,)’ ^3 0 



( 2 ) 


where ax, Uy, and az denote the Pauli matrices describing 
electron- and hole-like states {E/H). Likewise, ag = 
diag(5e,5h) is a 2 X 2-matrix in the space spanned by E 
and H and contains the effective (out-of-plane) g-factors 

and of the E and H bands, respectively. Like the 
other material parameters, the g-factors depend on the 
thickness d of the QW.— Whether the QW is in the 
normal or inverted regime, is determined by M. and B: 
li M./B <Q, the band structure is normal, whereas for a 
QW thickness d > dc, the band structure is inverted and 
M/B > 0. 

Moreover, we also investigate the effect of spin-orbit 
coupling (SOC) corrections, which—to lowest order—are 
described byii^ 


Hsoc 


^oreiA + ^ 

ri 


Tsiai 



+ BsiA2Py , 


( 3 ) 

where the first term describes the bulk inversion asymme¬ 
try of HgTe with its magnitude Aq, the remaining terms 
are the leading-order contribution to the structural inver¬ 
sion asymmetry due to the QW potential2ii2^ with the 
coefficient ^e, and the matrices are given by 


Tsiai 


0 icTp 
—icTp 0 


, rsiA2 


0 (Tp \ 

(Tp 0 y ’ 


Pbia 


0 -icTp \ 

icTp 0 y ’ 


( 4 ) 


= 0, spin is also a good quantum number and the 
spin-polarized eigenstates read as 


Kkir) 


fl(fLiy) 

y/L \ 9nkiy) 




( 7 ) 


where n is a band index, s the spin quantum number 
with its respective spinor Xs, L the length of the strip 
in the x-direction, and the functions and g^kiv) 

as well as the corresponding energy ens{k) can be deter¬ 
mined numerically or analytically for both cases (i) and 
(ii) from the respective Schrbdinger equations (see the 
Appendix lAl and Ref. 1^ for explicit solutions). If SOC is 
taken into account, the eigenstates are still translation- 
ally invariant, but no longer spin-polarized and thus are 
given by 


4'nfe(r) 



/ fLiy) \ 
yikiy) 
flkiy) 
\ 9 lkiy) / 


( 8 ) 


where n is again a band index, and we determine the 
functions /^^(y), 9nkiy)^ /nfc(2/)> and 172as well as 
the corresponding energy e„(fc) numerically with a finite- 
difference scheme. The eigenstates given by Eqs. dZD 
and (I5|) and their corresponding eigenenergies can then 
be used to calculate the magneto-optical conductivity via 
the Kubo formalism as will be discussed in the next sec¬ 
tion. 


B. Kubo formula for the magneto-optical 
conductivity 

Applying standard linear-response theory, one can 
write down Kubo formulas for the (magneto-)optical con¬ 
ductivities 


with (Tp = diag(l, 0). 

In this paper, we consider two geometries for the sys¬ 
tem described by Eqs. O-®: (i) bulk, that is, an in¬ 
finite system in the xy-plane, and—as the main focus 
of this work—(ii) a finite strip with the width w in the 
y-direction. For both cases, we apply periodic bound¬ 
ary conditions in the x-direction, and the confinement in 
case (ii) can be described by adding the infinite hard-wall 
potential 


CT;m(w) 


inL(^) 

Shu 


( 9 ) 


where the retarded current-current correlation function 
IV^{u) can be determined from the imaginary-time cor¬ 
relation function 


Him {iujn) = - / drfT /i(r)/m(0) 


( 10 ) 


0 for \y\ < w/2 
00 elsewhere. ' ' 

Thus, the total Hamiltonian reads as 

H = Ho + Hsoc + V{y)U (6) 

if SOC and confinement are taken into account. 

In both cases, (i) and (ii), translational invariance 
along the x-direction is preserved by the Hamiltonian ([5]) . 
Thus, the wave vector k in the x-direction is a good quan¬ 
tum number. Without SOC, that is, for Aq = 0 and 


via the formula H^^ (w) = Hj^ (w -|- i0+) and / and m 
denote the x- or y-directionsi^^— Here, Ii{t) denotes 
the charge current operator derived from the Hamilto¬ 
nian (O as described in the Appendix |B1 S the area of 
the QW, iun a bosonic frequency, r an imaginary time, 
T the imaginary time-ordering operator, (...) the thermal 
average, and /3 = l/{k^T) with the temperature T and 
the Boltzmann constant fee. 

Hence, we are left with the calculation of the retarded 
current-current correlation function, which can be deter¬ 
mined from Eq. (HOI) . In this paper, we investigate a sim¬ 
ple model: We assume that scattering by impurities can 








4 


be described by a constant, phenomenological scattering 
rate T/h and do not explicitly consider any other pro¬ 
cesses such as, for example, electron-phonon coupling.— 
Next, we introduce the spectral function, which in the 
spin-polarized case, that is, for Ag = 0 and = 0, is 
given by 


Ansik.U]) 


2hT 

[hio - ens{k) + -I- r2 ’ 


( 11 ) 


where e„s(fc) is the energy of the eigenstate labeled by 
the quantum numbers n, fc, and s of the system and p, is 
the chemical potential. 

If we insert the current operator in Eq. OT , express 
the Green’s functions in the resulting equation with the 
help of the spectral function d). calculate the sum 
over bosonic frequencies, and integrate over the resulting 
Dirac-(5 functions, we obtain the magneto-optical conduc¬ 
tivity tensors 


=737T 4n'(0,fc) 


A.'kSuj 


X Jduj'Ans{k,U}')An's{k,iO+Uj') 
X [riFDiflUj') — nFD(ftw' -I- huj)] 


( 12 ) 


and 


Im [cr;,j;(a;)] {^nn'(0: ^) Kn'(0> fc)]*} 

n,n\k,s 

X J dui' A ns {k,w')An's{k,Ul +uj') 

X [riFuikio') — n-pT>{hjjj' -I- hhj)] 

(13) 

where nFD(e) = 1/[exp(/3e)-I-1] and ctq = ^ jhw^ — 
Equations (IT^ and (fT^ also contain the dipole matrix el¬ 
ements (0, k) for transitions between bands n and n', 
where spin s and momentum k are conserved and which 
can be obtained from the current operator /; and the 
spin-polarized eigenstates in Eq. o as detailed in the 
Appendix IbI Finally, the conductivities Im[cri/(a;)] and 
Re[(7xy(w)] are then calculated from Eqs. m and m 
by using Kramers-Kronig relations. 

Equations (HB-dlSl) are written explicitly for the case, 
where spin is a good quantum number and optical tran¬ 
sitions are only permitted between states with the same 
spin quantum number. Spin conductivities af^{uj) can 
then also be defined by replacing the sum by ^ 
Eqs. (fT^ and (fOl) . In case SOC is considered and spin is 
no longer a good quantum number, the conductivities can 
be calculated in a similar way. Then, the equations de¬ 
termining the conductivities are given by Eqs. (HH)-®, 
but with the spin indices and the summation over spin 
omitted. Furthermore, the spin-unpolarized dipole ma¬ 
trix elements (0, k) for momentum-conserving transi¬ 
tions between bands n and n' obtained from the current 
operator Ii and the spin-unpolarized eigenstates ([S]) have 
to be used. 


In the following, we will use Eqs. (HU-CSl) and 
their spin-unpolarized generalizations to calculate the 
magneto-optical conductivities for geometries (i), that is, 
bulk, and (ii), that is, a finite strip. 


III. BULK 


As an introduction to the basic magneto-optical prop¬ 
erties of HgTe QWs, we first investigate the magneto¬ 
optical conductivity in a bulk HgTe QW without SOC, 
that is, where Ag = 0 and = 0 and spin s is a 

good quantum number. In this case, there is no con¬ 
fining potential V{y) in Eq. ([5]) and the eigenenergies are 
given by the dispersionless and highly degenerate LLs 
tns(k) = e„s, where n € Z is an integer. The LLs for 
n = 0 read as 

Cot = C + Ai - -p -h (14) 

and 

eoi = C-M - - yM bH, (15) 

while for n ^ 0 they read as 


2V\n\+B , 5e Affh c. , ^ ^ 

= C- - -^- - - ^bB + sgn(n) 




+ (m- 


p 


p 

'■B 


and 


(16) 


Cn4. — ^ 


2V\n\-B 

P 

'■B 


He 4“ Hh TJ i { \ 

—^— ^bB + sgn(n) 


2|n|A^ / 2B\n\-V 

II \ II 



2 


(17) 

for spin-up and -down electrons, respectively !^^ Here, 
n < 0 and n > 0 denote valence and conduction LLs, 
respectively, while of the two LLs at n = 0, one belongs to 
the conduction band and the other to the valence band.— 
The critical magnetic field, where those two zero levels 
cross and which separates the QSH and QH regimes, can 
be calculated from the condition eg-|- = egj, and yields^ 


Be 


_ M _ 

eBjh - (5e +5h)HB/4' 


(18) 


Calculating the dipole matrix elements d^^/(0, fc) from 
the corresponding eigenstates yields the optical selection 
rules 


|n| —>■ |n ± 1|, s —!► s, fc —!► fc 
for transitions between LLs. 


(19) 
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FIG. 3. (Color online) Real and imaginary parts of the (a) lon¬ 
gitudinal magneto-optical and (b) optical Hall conductivities 
Oxy(uj) in a bulk HgTe QW of thickness d = 7.0 
nm with some transitions explicitly labeled. 

To illustrate the resulting absorption spectrum, Fig. [3] 
shows the real and imaginary parts of the numerically^^ 
obtained magneto-optical conductivities crxxi<-^) and 
axy(uj) for the material parameters A = 364.5 meV nm, 
B = —686.0 meV nm^, C = 0, V = —512.0 meV nm^, 
M = —10.0 meV, = 22.7, and gi, = —1.21, correspond¬ 
ing to a QW thickness d = 7.0 nm > dcr^ that is, for 
parameters in the QSH regime (at B = 0). The magnetic 
field is chosen to be B = 5 T, that is, a magnetic held 
below Be ~ 7.4 T for which the model still shows an in¬ 
verted band structure (see Fig. [I, while the remaining 
parameters are chosen to be T = 1 K, ^ = 8 meV, and 
F = 1 meV. Moreover, we note that the remaining com¬ 
ponents of the conductivity tensor can be determined 
from ayy{uj) = axx(<^) and axy(oj) = —cTyx(uj) for the 
bulk system considered in this section. 

As can be seen in Figs. [3] (a) and (b), there are multi¬ 
ple peaks in the absorptive components Re [<yxx{^)\ and 
luY[axy{ui)\ corresponding to transitions between an oc¬ 
cupied and an unoccupied LL state where the selection 
rules CH) have to be satished. Here, both intraband tran¬ 
sitions, that is, transitions between only conduction LLs 
or only valence LLs, at low energies and interband transi¬ 
tions, that is, transitions between valence and conduction 
LLs, at higher energies are possible. 

Due to the non-linear dependence of the LLs m- 
(II3 on n, different allowed transitions between LLs have 
different energies resulting in the multiple absorption 
peaks shown in Fig. [3] (a)i^ This behavior, also ob¬ 
served experimentally^ and reminiscent of the situation 
in graphene^^— differs markedly from the behavior in 
a normal 2D electron gas, where equidistant LLs lead to 
only one absorption peak. Moreover, we find that, while 
every transition satisfying Eq. (US occurs, different tran¬ 
sitions are most probable for spin-up and spin-down LLs: 
If n > 1, the squared dipole matrix element for a tran¬ 
sition —n —^ n -|- 1 is typically five to ten times larger 
than the one of — n —> n — 1 for spin-up LLs and vice 



FIG. 4. (Golor online) Real part of the longitudinal magneto¬ 
optical conductivity in a bulk HgTe QW of thickness 

d = 7.0 nm for different (a) temperatures T, (b) chemical 
potentials g, and (c) magnetic fields B. 


versa for spin-down LLs. Consequently, the prominent 
interband transition peaks are governed by —n n + 1 
for spin-up LLs and by —n —>■ n — 1 for spin-down LLs 
in Figs. [3] (a) and (b), while other interband transitions 
contribute much less to the absorption spectrum. 

In addition to the absorptive components, the refrac¬ 
tive components Im [cr^xC^)] and Re [axy{uj)] are also dis¬ 
played for completeness. Both, the absorptive and refrac¬ 
tive components, are generally needed to calculate optical 
observables, such as reflection coefficients and angles, for 
example. 

In Fig. m the dependence of Re [crxx(a;)] as a function 
of the frequency on different parameters is displayed for 
a bulk HgTe QW again in the QSH regime (at B = 0). 
The temperature dependence is illustrated in Fig. |4] (a), 
which shows Re [cTxx{<^)] for a fixed magnetic field, chemi¬ 
cal potential, and broadening. Here, the main feature ob¬ 
served is that with increasing T, additional (intraband) 
transitions (forbidden at T = 0 ) become more proba¬ 
ble which then give rise to new peaks—mainly at low 
frequencies. For higher frequencies, on the other hand, 
the magneto-optical conductivity remains largely unaf¬ 
fected, although the interband peaks are slightly reduced 
as spectral weight is transferred to lower energies, while 
the total spectral weight is conserved. 

Figure m (b) shows Re[crxx(w)] for several different 
chemical potentials and a fixed magnetic field, temper¬ 
ature, and broadening. As g increases, the intraband 
transition peaks move to lower energies, while the gap 
between intraband and interband transitions increases. 
Likewise, with decreasing B, the energy of the intraband 
transition peaks decreases and tends to zero, as can be 
seen in Fig. S] (c), which displays Re [(Jxxi^)] for several 
different values of B and fixed g, T, and F. 
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FIG. 5. (Color online) Real parts of the longitudinal and 
spin longitudinal magneto-optical conductivities, and 

(cu), respectively, in a bulk HgTe QW of thickness d = 7.0 
nm (a) at the critical magnetic field B = ~ 7.4 T and (b) 

at i? = 5 T. In panel (a), /r is chosen to be at the energy 
of the degenerate spin-up and spin-down zero mode Landau 
levels, while in panel (b) ^ is chosen to be at the energy of 
the spin-up zero mode Landau levels. 


The behavior of the intraband transition peaks with 
decreasing magnetic field or with increasing chemical po¬ 
tential can be qualitatively explained as originating from 
the LL spectrum in the vicinity of the Fermi level: For a 
fixed magnetic field, the LL spacing near the Fermi level 
decreases if the absolute value of the chemical potential 
is increased and thus situated in a denser region of the 
LL spectrum. Consequently, the energies of the intra¬ 
band transitions are also decreased. On the other hand, 
with decreasing magnetic field the LL spacing decreases 
giving rise to lower energies of the intraband transitions. 
Moreover, the amplitudes of the interband peaks decrease 
with decreasing magnetic fields, as seen in Fig.|3](c). This 
can be interpreted as arising from optical transitions be¬ 
tween increasingly denser regions of the LL spectrum, 
where the energy difference between different transitions 
is small compared to the broadening due to scattering. 

In the above discussion, we have investigated a QW 
with parameters corresponding to the topological regime 
(at B = 0). However, the above conclusions on the be¬ 
havior of the magneto-optical conductivity also apply to 
QWs with a thickness d < dc, that is, QWs in the topo¬ 
logically trivial regime. We conclude our discussion of 
bulk HgTe QWs by investigating also the spin magneto¬ 
optical conductivity (T^^®(a;), where one can indeed find a 
signature of the transition between the QSH and the QH 
regimes: As shown in Fig. [21 the transition between the 
QSH regime at finite magnetic field and the QH regime 
in QWs with a thickness d > dc occurs when the spin-up 
and spin-down zero LLs cross at a critical magnetic field 
Be- 

If the sample is undoped as chosen in Fig. [5] (a), the 
chemical potential lies between the two zero LLs. This, 
however, means that for B = Be the chemical potential 


will lie in the two degenerate zero LLs. Thus, there is 
an additional Drude transition peak at zero frequency 
as shown in Fig. |S] (a). This peak originates from two 
different transitions, namely spin-conserving transitions 
within the spin-up or spin-down zero LLs. As the squared 
dipole matrix elements for transitions within the spin-up 
and spin-down zero LLs are the same, both transitions 
contribute equally to the absorption peak at Sw = 0. 
However, this means that there is no peak at hu) = 0 for 
the spin conductivity as spin-up and spin-down 

transitions cancel each other exactly. 

This only happens at the critical field Be in an undoped 
sample and is different from the situation if /r in a doped 
sample crosses one LL, an example of which is shown in 
Fig. |5| (b) for comparison. In this case, there is also a 
Drude peak at fiw = 0 for axx{^), but since this peak 
arises from transitions within only one LL with a given 
spin quantum number, (oj) also exhibits a peak at 
zero frequency and can thus be distinguished from the 
peak arising from the two zero LLs at Be- 


IV. FINITE STRIP 

While in the previous section we have investigated the 
magneto-optical properties of bulk HgTe QWs without 
SOC, which did not account for the presence of edge 
states and where we did not find significant differences 
between the QSH and the QH regimes, we will now turn 
to a finite-strip geometry described by the confining po¬ 
tential in Eq. ©• Since the finite-strip geometry contains 
boundaries, both chiral QH edge states as well as helical 
QSH edge states emerge at these boundaries under the 
appropriate conditions: If the width w of the finite wire 
is large compared to Zb, well localized QH edge states 
form in addition to the bulk LLs. In the QSH regime, 
that is, in a QW with d > dc and below Be, there are 
also QSH edge states in addition to the QH states. 

This can also be seen in Fig. |6l which shows the real 
parts of the magneto-optical conductivities o'xx(w) and 
ayy{uj) for a finite-strip of width w = 200 nm for differ¬ 
ent QW thicknesses and magnetic fields corresponding to 
different regimes^i without SOC (Aq = 0 and = 0): a 
QW with d < dc, that is, a system with no QSH states 
at any magnetic field [Fig. [5] (a)], a QW with d > dc and 
B < Be, that is, a system possessing QSH states at the 
edges of the strip [Fig.|n](b)], and a QW with d > dc, but 
B > Be, that is, a system with no QSH states [Fig.|n](c)]. 
For each of the regimes, we have chosen /r to lie at the 
neutrality point between the uppermost valence LL and 
the lowest conduction LL. The insets of Figs. [HI (a)-(c) 
show the respective energy spectra and the positions of 

In the regimes without any QSH edge states [Figs.|n](a) 
and (c)], the band structure is gapped as every electron¬ 
like LL is a conduction band and situated above the hole¬ 
like valence bands. The edge states associated with an 
electron-like LL have positive curvature, while those as- 







7 



FIG. 6. (Color online) Real parts of the longitudinal magneto¬ 
optical conductivities (Jxxioj) and ayy{u)) in (a) a HgTe QW 
of thickness d = 5.5 nm and with a magnetic field B = 5 T, 

(b) a HgTe QW with d = 7.0 nm and B = 5 T, and (c) a 

HgTe QW with d — 7.0 nm and B = 10 T. Here, spin-orbit 
coupling corrections are not included. The solid orange lines 
represent = (Tyy{(jj) in a bulk system, while the solid 

black and the dashed green lines represent (Jxx{oj) and (Jyy{uj) 
in a finite strip of width w = 200 nm. The insets in panels (a)- 

(c) illustrate the respective energy spectra of the finite strip 
and the positions of the chemical potential (dotted lines). 


sociated with a hole-like LL have negative curvaturei^ 
Figure m (b), on the other hand, shows a different situa¬ 
tion as the lowest electron-like (spin-up) LL is a valence 
band and the highest hole-like (spin-down) LL is a con¬ 
duction band. Thus, there is a crossover between the 
dispersions of electron- and hole-like bands and one con¬ 
sequently finds counterpropagating, spin-polarized QSH 
states, in addition to QH edge states propagating in the 
same direction at a given boundary regardless of spin. 
For comparison, the magneto-optical conductivities for a 
bulk HgTe QW as investigated in Sec. Hill are also dis¬ 
played in Fig. ini for the same magnetic fields and QW 
thicknesses d as the finite-strip structures. 

Compared to the situation in a bulk HgTe QW, we 
can see that in a finite strip there are still pronounced 
absorption peaks arising from the bulk LLs. However, the 
spectral weight of these peaks is reduced as transitions 
between edge states are also possible. Since the energy 
of these transitions can differ from the energy difference 
between two bulk LLs and the total spectral weight is 
conserved, some spectral weight is transferred from the 
peaks to the regions between the absorption peaks, a 
phenomenon clearly seen in Fig. [51 

A second difference compared to the bulk system is 
that, even if the chemical potential is not situated in a 
bulk LL, there is a Drudc-like absorption peak at low en¬ 
ergies for cjxx{^)- Analogous to the situation in Fig. [SI 
where a Drude peak at low energies originated from tran¬ 


FIG. 7. (Golor online) Real part of the total longitudinal 
magneto-optical conductivity of a finite strip with 

width w = 200 nm, F = 1 meV, and chemical potential fj. — 8 
meV at temperature T = 1 K as a function of the inverse 
magnetic field 1/B and fkd for band parameters correspond¬ 
ing to a HgTe QW of thickness d = 7.0 nm without spin-orbit 
coupling corrections. The vertical line indicates the inverse 
critical magnetic field separating the QSH and QH regimes 
and as calculated by Eq. m- 


sitions within bulk LL states though, this absorption 
peak arises from transitions that occur when the energy 
of a QH or QSH edge state is at the chemical potential 
or close to it (within the broadening). This mechanism is 
present in both the QH and QSH regimes if the chemical 
potential is above or below the bulk band gap. The dif¬ 
ference between those two regimes, however, is that edge 
states exist at every energy in the QSH regime, while 
there are no states with energies inside the bulk gap in 
the QH regime (see the insets in Fig. 0). Hence, if in the 
QH regime the chemical potential lies inside the gap and 
the gap exceeds the energy scale associated with broad¬ 
ening, no low-energy absorption peak occurs as seen in 
Figs. [ 5 ] (a) and (c). Moreover, we note that this Drude 
peak appears only along the unconfined direction, but 
not for ayy{ui) along the confined direction. This is to be 
expected because only the x-direction is associated with 
free acceleration. 

The disappearance of the Drude peak if ^ remains in¬ 
side the bulk gap during the transition from the QSH 
to the QH regime is also illustrated in Fig. [T] Here, the 
absorption spectrum of a finite-strip geometry for a QW 
with d = 7.0 nm > dc (and without SOC) is shown as 
a function oi 1/B and Huj at low temperatures and fixed 
chemical potential. For 1/B = 0.2/T, the situation is de¬ 
scribed by Fig.[n](b). As 1/B is decreased, that is, as the 
magnetic field is increased, the absolute value of the LL 
energies and thus the photon energies of their associated 
absorption peaks increase. At low photon energies, there 
is a finite absorption, which vanishes for magnetic fields 































































FIG. 8. (Color online) Real parts of the longitudinal magneto¬ 
optical conductivities a^xioj) and ayy{uj) of a finite strip with 
width w = 200 nm, a magnetic field B = 5 T, and band 
parameters corresponding to a HgTe QW of thickness d = 7.0 
nm with and without spin-orbit coupling corrections as given 
in Eq. ©■ The respective energy spectra and the position of 
the chemical potential (dotted line) are shown in the inset. 


above the critical field Be when the transition from the 
QSH into the QH regime occurs. 

We now address the effect of SOC. While the states 
in Fig. [6] have been perfectly spin-polarized and char¬ 
acterized by the spin quantum number s, the situation 
changes if SOC is considered, that is, if Aq and/or are 
finite in the total Hamiltonian ([3]). This is illustrated in 
Fig.l where the real parts of the magneto-optical con¬ 
ductivities Uxxi^) and ayy(u)) are displayed for a QW 
of thickness d = 7.0 nm > dc at H = 5 T (see above) 
with the additional SOC parameters Aq = 1.6 meV and 

= 16.0 meV nmiii The inset of Fig.[8]compares the en¬ 
ergy spectrum near the band gap in the presence of SOC 
with the spin-polarized energy spectrum in the absence 
of SOC. The most striking effect appears at the cross¬ 
ing between the spin-up and spin-down states, where a 
SOC gap is opened up. Consequently, the low-energy ab¬ 
sorption peak present in the QSH regime without SOC 
is reduced if the chemical potential is situated inside this 
small gap. However, as long as the gap opened by SOC 
does not significantly exceed the broadening F, an in¬ 
creased absorption for c7xx(<^) can still be observed in the 
QSH regime at low photon energies hut. As bands farther 
away from the zero LLs are affected even less by SOC 
for the parameters given above, the absorption spectrum 
at higher energies remains nearly unaltered compared to 
spin-polarized case. 

Next, we turn our attention to the effect of tempera¬ 
ture on the absorption spectrum in general and on the 
Drude-like peak arising from QSH states in particular 
and compare this to the QH regime. For this purpose. 
Figs, ini (a) and (b) show Re [cTxxi^)] in a finite-strip ge¬ 
ometry (without SOC) at different temperatures for a 
QW with d = 7.0 nm > dc in the QSH and QH regimes, 



FIG. 9. (Golor online) Real part of the longitudinal magneto¬ 
optical conductivity of a finite strip with width w = 

200 nm, band parameters corresponding to a HgTe QW of 
thickness d — 7.0 nm without spin-orbit coupling corrections, 
and magnetic fields of (a) B = 5 T and (b) B = 10 T for 
different temperatures T. The energy spectra and the posi¬ 
tions of the chemical potential (dotted lines) are shown in the 
insets. 



T[K] 


FIG. 10. (Color online) Real part of the total longitudinal 
magneto-optical conductivity of a finite strip with 

width w = 200 nm as a function of the temperature T and hw 
for B = 5 T, ^ = 8 meV, F = 1 meV, and band parameters 
corresponding to a HgTe QW of thickness d = 7.0 nm without 
spin-orbit coupling corrections. 


respectively, with the chemical potential situated be¬ 
tween the zero LLs. Complementary, Fig. (TU] displays 
Re [(Jxx{^^)\ for the setup of Fig. [3] (a) as a function of T 
and huj. 

Similar to the situation in a bulk QW illustrated in 
Fig. 0](a), additional transitions especially, but not ex¬ 
clusively, at low energies become more probable. This 
results in additional absorption peaks and an increase 
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FIG. 11. (Color online) Real parts of the (a) total and (b) 
spin longitudinal magneto-optical conductivities aj:x(oj) and 
of a finite strip with width w = 200 nm at low fre¬ 
quencies ujo as a function of the temperature T for different 
setups. 

of the absorption at low energies, while spectral weight 
is transferred away from the intraband peaks as can be 
seen in Figs. and ITOl Since the spin-up and spin-down 
zero modes are very close to /j- in these setups, additional 
transitions involving one of these LLs become especially 
likely if T is increased and consequently the LL below /i 
is depopulated, while the LL above /i is populated. This 
is particularly striking in Fig. IH] (b), where already at 
T = 100 K pronounced peaks appear due to transitions 
n = 0 —>■ 1 and — 1 —>■ 0 for spin-up and spin-down LLs, 
respectively. 

As the absorption at low energies becomes more prob¬ 
able, an enhancement of the low-energy absorption peak 
in the QSH regime can be observed in Figs. [S] (a) and [TUI 
On the other hand. Fig. (b) shows that such a peak 
also appears in the QH regime with increasing temper¬ 
ature. The smaller the gap between conduction and va¬ 
lence LLs, the smaller is the temperature for the onset 
of this effect. Thus, at high T both regimes exhibit a 
low-energy absorption peak for axx{<^)- 

This is also corroborated by Fig.[TT](a), where the real 
part of at low energies is displayed as a function 

of T for a QW in the QSH regime, in the QSH regime 
with SOC, and in the QH regime with the chemical po¬ 
tential situated between the conduction and valence LLs. 
At low temperatures, there is a finite absorption in the 
QSH regime, even if reduced by SOC, whereas there is 
no absorption in the QH regime. With increasing T, 
the absorption also increases, in both the QSH and QH 
regimes, although the increase is more pronounced in the 
QSH regime. Moreover, SOC corrections are less pro¬ 
nounced at higher temperatures. Without SOC, a spin 
conductivity tr^jf (w) can be defined, whose temperature 
dependence is shown in Fig. [TT] (b) and closely follows 
the temperature dependence of axx{<^)- 

Finally, we study the dependence of the low-energy 
absorption on the chemical potential. For this purpose. 



T ^ ^ ^ ^ ^ ^ I 

0 5 10 15 20 


(1 [meV] 

FIG. 12. (Color online) Real parts of the (a) total and (b) 
spin longitudinal magneto-optical conductivities and 

cr^)f(cj), respectively, of a finite strip with width w = 200 nm 
at low frequencies cjo as a function of the chemical potential 
p for different setups. 

Fig. HU] (a) shows the real parts of the conductivity (Jxx{<jj) 
at low hu! and T as functions of /i for a QW in the QSH 
regime, in the QSH regime with SOC, and in the QH 
regime. The corresponding spin conductivities cr^jf (w) 
for the regimes without SOC are displayed in Fig.lTUKb). 
As can be seen in Fig.[T2|(a), the absorption in the QSH 
regime is finite and usually higher than in the QH regime. 
Spin-orbit coupling in the QSH regime only has a signifi¬ 
cant effect close to the crossing between the electron- and 
hole-like edge dispersions, that is, in the energy interval 
approximately between 5 meV and 15 meV. In this re¬ 
gion, o-xx(w) is reduced due to SOC, although it does not 
vanish completely. The behavior with varying chemical 
potential in the QH regime, on the other hand, is differ¬ 
ent: If fi is situated in the gap between the valence and 
conduction LLs, there is no Drude-like peak and axx{^) 
vanishes. 

Thus, by sweeping the gate-voltage and thereby vary¬ 
ing the position of the chemieal potential, it is possible 
to distinguish between the QSH and QH regimes on the 
basis of their their low-energy absorption as depicted in 
Fig. [Ill (a). We note that this behavior of the Drude 
peak in the absorption spectrum is consistent with and 
reflects the expected behavior of the dc conductance in 
the QSH and QH regimes^ 

This is also substantiated in Fig. |T3| which shows the 
low-energy absorption of a QW of thickness d = 7.0 nm 
without SOC corrections as a function of both the chem¬ 
ical potential as well as the magnetic field. While in the 
QSH regime, B < Be ~ 7.4 T, Re [cjxxi^)] never vanishes 
entirely, there is a region where Re {uxxij^)] is completely 
suppressed in the QH regime above B^. With increasing 
magnetic field, the extent of this region grows as the gap 
between conduction and valence-LLs increases. Outside 
the region of suppressed conductance and low-energy ab¬ 
sorption, the spectrum displays an oscillatory behavior 
as seen in Figs. [TUI and [TUI 
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FIG. 13. (Color online) Real part of the total longitudinal 
magneto-optical conductivity axa:{io) of a finite strip with 
width w = 200 nm and F = 1 meV at low energies fiojo = 0.1 
meV and temperature T = 1 K as a function of the chemical 
potential jj, and magnetic field B for band parameters cor¬ 
responding to a HgTe QW of thickness d = 7.0 nm without 
spin-orbit coupling corrections. The vertical line indicates the 
critical magnetic field separating the QSH and QH regimes 
and as calculated by Eq. m- 


the QSH and QH regimes based on these LL peaks. 

We find that these bulk LL peaks also occur in the 
finite strip, where the presence of edge states results in 
an additional Drude-like absorption peak in cr^x orig¬ 
inating from low-energy transitions at the Fermi level. 
This Drude peak is always present in the QSH regime, 
while it vanishes in the QH regime if the chemical poten¬ 
tial is situated in the bulk gap. If SOC corrections are 
included, we find that their effect is to open up a small 
gap in the QSH spectrum, which leads to a reduction, 
but not to a complete disappearance of this Drude peak. 
By sweeping the gate-voltage of a finite QW structure 
and thereby varying the position of the chemical poten¬ 
tial, it is therefore possible to experimentally distinguish 
between the QSH and QH regimes on the basis of the sta¬ 
bility (QSH) or disappearance (QH) of the Drude peak. 
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Einally, we remark that our description does not 
take into account many-body effects, such as (edge) 
magneto-plasmon resonances which can potentially play 
a prominent role for the absorption especially in narrow 
strips4^“— Since the scheme that we propose to distin¬ 
guish between the QSH and QH regimes is based on the 
dc conductivity, however, we think this scheme to be ro¬ 
bust, even in the presence of edge magneto-plasmons, al¬ 
though the absorption at low, but finite energies can ob¬ 
tain a richer structure than the single-particle absorption 
calculated in this paper. Moreover, our single-particle 
calculations provide a useful benchmark against which 
to test experimental results. 

V. CONCLUSIONS 


Appendix A: Eigenspectrum and eigenstates 
1. No spin-orbit coupling corrections 

Inserting the ansatz from Eq. 0 in the Schrbdinger 
equation given by the Hamiltonian (jS]) for Ag = Ce = 0 
and introducing the transformation 

C = ay) = V2{y-llk)/l, (Al) 

yields a system of two differential equations to determine 
the eigenenergy E = ens{k) and eigenstates with given 
spin and momentum quantum numbers s and k. If the 
y-components in Eq. © are defined as f^^{y) = /«[((?/)] 
and glk^y) = 5s[C(2/)], the system is given by 


Inverted HgTe/CdTe QWs exhibit helical QSH states 
below a critical magnetic field B^, above which the band 
order is normal and only QH states can be found. In this 
work, we have studied the magneto-optical properties 
of such HgTe/CdTe QWs by calculating the magneto¬ 
optical conductivity using linear-response theory. We 
have considered both an infinite 2D system as well as 
a finite strip, that is, a 2D system that is confined in one 
direction. The normal and the inverted regimes both ex¬ 
hibit a series of pronounced absorption peaks correspond¬ 
ing to different transitions, both intraband as well as in¬ 
terband transitions, between non-equidistant LLs in the 
infinite system. Thus, it is hard to distinguish between 


0 = 

+ 



St +9c) /s(C) 


Next, we impose boundary conditions 
direction, namely 


( gJsiC) \ 

V 5hffs(C) J ' 

(A2) 

along the y- 


lim fsiO= lim g^((/) = 0 (A3) 

(^—>•±00 
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for the infinite bulk geometry (i) and 

fs [C(±«^/2)] = ~9s [C(±«^/2)] = 0 (A4) 

for the finite-strip geometry (ii). 

We use a finite-difference scheme to solve the system 
of differential Eqs. (jA2() numerically. In addition, we also 
apply a second method to obtain and check the solutions 
of Eq. (IA2I) : As detailed in Ref. |4^ and similar to the 
procedure in Ref. [s^ the general solution of Eq. (ESI) 
can be written down in the form of parabolic cylindrical 
functions. By invoking the appropriate boundary con¬ 
ditions on this general solution, one can then obtain a 
transcendental equation to determine the eigenspectrum 
and eigenstates. The boundary conditions for a bulk sys¬ 
tem given by Eq. (IA3|) . for example, lead to the LLs in 
Eqs. (O-dlZl) and their corresponding eigenstates, where 
the parabolic cylindrical functions are reduced to Her- 
mite polynomials.— 


2. Spin-orbit coupling corrections 


If SOC corrections Aq and/or are considered and 
the ansatz from Eq. m as well as the transforma¬ 
tion in Eq. (EH) are used, we redefine the functions 

fnkiy) = /i[C(y)], 5nfc(y) = MCiv)], flkiv) = /2[C(y)], 

and glk(y) = MCiy)]- Then, the eigenenergy E = e„(fc) 
of an eigenstate with momentum k is determined by the 
system of four differential equations 



f Mo \ 

ffi(C) 

Mo 


\M 0 J 


/ Mo \ 

-hiO 

Mo 

\ -MO / 


+ 




02A 

Ib 


^ U-dcVgiiO ^ 
ll + dc)Mo 
“ (§ + MO 

V - (2 -^c) Mo ! 


I gJiiO \ 

y-sB g^MO 

2 -g^MO 
V -gi^MO / 


( -MO \ 


( {l+dc)M 0 \ 

Mo 

iv^Ce 

0 

5 i (0 

Ib 

'"A 

1 

1 

V -Mo y 


V o' / 


(AS) 

and the boundary conditions Eqs. (IA3I) or Eqs. (IA4I) . 
where s is to be replaced by the indices 1 and 2. The 
system of differential Eqs. (lASp is then solved numerically 
with a finite-difference scheme. 


Appendix B: Current operator and dipole matrix 
elements 


To calculate the optical conductivity via the Kubo for¬ 
mula given by Eqs. m and (113, we need to know the 
charge current operator I. In the presence of an arbi¬ 
trary magnetic vector potential A(r), Eqs. ([T|) and ([3]) 
describing the HgTe QW have to be generalized by 


Hq —CI 4 -l- — 


VIa 


BY 5 




Mx + ^l) 


ATi 


AY 


2 . [V X A(r)] Tg 

H--- 


(Bl) 


and 


Hsoc = AoTbia + (Bsiai^x + rsiA2d‘y), (B2) 

respectively. Here, the kinetic momentum operator tt = 
p -I- eA(r), Tg = (rg,r|,rg) and the additional 4x4 
matrices 


-pX _ 

^ 9 ~ 


0 5l|l2 
5||12 0 


j^y — 
’ 9 


0 -igi|l 2 

igiil 2 0 


(B3) 


as well as the effective in-plane g-factor have been 
introduced— 

For an arbitrary (normalized) state 4'(r), the corre¬ 
sponding energy expectation value of the total Hamilto¬ 
nian iJ, given by Eq. (ED and generalized by Eqs. m 
and (IB2I) . as a functional of the vector potential A(r) can 
be obtained as 

E[A]=Y,J (B4) 

where the sums over a and /3 refer to the four bands con¬ 
sidered, that is, \H t), \E i), \E[ |). The particle 

current density j(r) of this state 4'(r) can be determined 
by a variational method: 

5E = E[A + 5A]- E[A] = e j dV j(r)(5A(r). (B5) 

This procedure yields the probability current density 
j(r) = je(r) +ji(r) composed of the external current den¬ 
sity, 

je(i-) = 




EU 2?(i4)„^ + e(r5)„^ [4/* (Vd/^) - (vvk*) vh^] 

a/3 I L J 

+ x (r)aa (l4)„a + B (Ts) 


3a/3 


-T (BsIa)„; 3 4'a4'/3 


(B6) 

where T = (ri,r2,0) and Tsia = (Tsiai, rsiA2, 0 ), and 
the internal current density. 


J.W = X 


aj3 


(B7) 
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As we are dealing with a 2D system, Eqs. (IB6I) and (lB7l) 
are to be read as applying only to the x- and y- 
components. 

We note that the external current density given by 
Eq. (IB6|) could also have been obtained by calculating 


the velocity operator v = r,H /ih and using j = 

(v^*) + (v'l')* 'h] /2. Since for the parameters and 
magnetic fields investigated in this work, the internal cur¬ 
rent is very small compared to the external current, we 
neglect its contribution and only consider the external 
current density from now on. 

By promoting the wave functions 'l'a(r) and 4'/3(r) in 


Eq. (IB6p to field operators 'l'a(r) and 'I',a(r), we obtain 
the current density operator je(r). As basis for the field 
operators, we use the eigenstates of the respective system, 
that is, the states given by Eq. 0 and determined from 
Eq. dAH) for a system where spin is a good quantum 
number and the states given by Eqs. dH]) and (IA5I) if SOC 
corrections are taken into account. The charge current 
operator I is then obtained from 


i = 



which yields the components 


k,s,n,n' 


(B8) 


(B9) 


il = -e Y (BIO) 

k,n,n' 


in the /-direction for the case with and without SOC cor¬ 
rections. Here, the operators {cnks) and (c„fc) 
create (destroy) an electron in a state given by Eqs. ([7]) 
and dll), respectively. The corresponding dipole matrix 
element d„*/(0,/c) is calculated by replacing 4'^(r) and 
'l'^(r) in Eq. (IB6I) with the eigenstates [^'^fc('^)]a 
[^n'fc (*^)]/3 from Eq. dZ]) and integrating over d^r. In the 
case of an infinite bulk system, analytical formulas can be 
derived for the dipole matrix elements using the eigen¬ 
states given in Ref. However, those dipole matrix 
elements are quite cumbersome and not particularly elu¬ 
cidating. If SOC corrections are considered, the dipole 
matrix element (0, k) is determined in a similar way 
from the eigenstates [ll'nfc(i‘)]a and [\I'n/fc(r)]o of Eq. dll)- 
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